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Maximum-likelihood estimation is applied to identification 
of an unknown quantum mechanical process represented by 
a "black box". In contrast to linear reconstruction schemes 
the proposed approach always yields physically sensible re- 
sults. Its feasibility is demonstrated using the Monte Carlo 
simulations for the two-level system (single qubit). 

PACS number(s): 03.65. Bz, 03.67.-a 



During recent years great attention has been devoted 
to the measurement of quantum, state of various simple 
quantum mechanical systems. All proposed reconstruc- 
tion techniques follow the common underlying strategy: 
A set of measurements is performed on many identically 
prepared copies of the quantum state which is then es- 
timated from the collected data. Feasible reconstruction 
schemes were devised for a wide variety of systems includ- 
ing the modes of running electromagnetic field (optical 
homodyne tomography |l|,Q] and unbalanced homodyn- 
ing p]), cavity electromagnetic field Hg, motional state 
of ion in Paul trap [g[0] , vibrational state of the molecule 
lH and spin Q. 

These significant achievements stimulated develop- 
ment of a new remarkable branch of the reconstruction 
techniques that allow for the experimental determination 
of the unknown quantum mechanical processes (M-nSl. 
This is of great practical importance, because such a tech- 
nique may be used, e.g. to evaluate experimentally the 
performance of the two-bit quantum gate - a building 
block of quantum computers [|0| . The usual set-up con- 
sidered also in this paper is shown in Fig. 1. Note that 
a similar experimental configuration can also allow for a 
complete characterization of quantum measurement p6[ . 
The input state prepared by an experimentalist and char- 
acterized by a density matrix gin enters the "black box" 
where it is transformed into the output gout • The task for 
the experimentalist is to retrieve information on the phys- 
ical process hidden in the black box from the measure- 
ments on the output states Pout- The only assumption 
taken for granted here is that the mapping pout = GQm is 
linear, as dictated by the linearity of quantum mechanics. 
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Here gtj — {i\g\j) are density matrix elements in some 
complete orthogonal basis of states spanning the Hilbert 
space on which the density operator g acts. As illus- 
trated in Fig. 1, the system may be entangled with the 
environment and the transformation G need not preserve 
purity of the state. The Green superoperator Q can de- 



scribe a diverse variety of the physical processes, such as 
unitary evolution, damping and decoherence. From the 
reconstructed superoperator Q one may further estimate 
the Liouville superoperator C, which governs the evolu- 
tion of density matrix in the black box, g = Cg. If the 
superoperator C exists, then G = exp(£T), where r is the 
interaction time, and an inversion of this relation yields 

The estimation of the elements Q^j by means of linear 



algorithms has been addressed in several papers [|10|-[12[ . 
Provided that mapping between known input and output 
states is given, the unknown parameters Q^j may be ob- 
tained from (m using the system of linear equations. This 
linear reconstruction procedure is simple and straightfor- 
ward, but it suffers from one significant drawback. The 
elements 5f' are estimated as a set of seemingly unrelated 
numbers. However, Q^j cannot be arbitrary because the 
linear mapping Q must preserve the positive semidefinite- 
ness and trace of the density matrix. These conditions 
impose bounds on the allowed values of Gfj ■ In this Rapid 
Communication the superoperator G is reconstructed us- 
ing maximum-likelihood (Max-Lik) . It allows to incorpo- 
rate naturally all the constrains of quantum theory. Since 
one can only collect a finite amount of data, the linear 
mapping cannot be determined exactly. In accordance 
with the probabilistic interpretation of the quantum the- 
ory, the Max-Lik estimation answers the question "Which 
process is most likely to yield the measured data?". How- 
ever, Max-Lik solution is not only an estimation, but 
represents a genuine quantum measurement associated 
with a quantum device. 

Black box 




FIG. 1. Sketch of experimental set-up for determination 
of the quantum-mechanical process. The input state £iin is 
prepared in the preparator P and enters the black box where 
it is transformed to the output state Qout = GOin which may 
be entangled with the environment. The detector D measures 
some observable of the output Qout- 



Due to its nonlinearity, the Max-Lik estimation is com- 
putationally much more expensive task then the linear 
procedures. This is the prize paid for the physically 
sound result. Max-Lik estimation has been applied to 
various problems recently: To the measurements of the 
quantum phase shift |1^, a coupling constant between 
atom and a cavity electromagnetic field |iaL or the pa- 
rameters of quantum-optical Hamiltonian|19| . Recon- 
struction of generic quantum state using the Max-Lik 
estimation and its interpretation as quantum measure- 
ment has been proposed in [Eoj . Subsequent Monte Carlo 
simulations, performed for the quantum states of electro- 
magnetic field modes and spin pl|-p3[|, illustrated a fea- 
sibility of this technique. Here we shall demonstrate that 
the Max-Lik estimation is also suitable for determination 
of the generic quantum mechanical processes. 

The sought superoperator Q can be found as that one 
maximizing the likelihood function C[G]. Let us consider 
n measurements described by positive operator-valued 
measures (POVM) H^™), m = 1, . . . ,n. Then C[g] reads 
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Thus X is positive semidefinite hermitian matrix. This is 
the desired condition revealing a domain of the allowed 
parameters Q^^ (or, alternatively, Xij)- The matrix x is 
parameterized by N'^ real numbers, but the condition (3) 
imposes N"^ real constraints so that the number of inde- 
pendent parameters reads N'^ — N"^ . When the operator 
expansion (H) is substituted into Eq. (0), one obtains 
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where a^^„ — {m\A^^Aj\n). The constraints can be also 
expressed in terms of Q^!- , 
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where fm is (relative) frequency for detection of 11^™^ 
Maximum of this function should be found in the do- 
main of physically allowed superoperators G, whose de- 
termination is crucial for successful implementation of 
the Max-Lik estimation. The linear positive map (Q) can 
be conveniently cast into the form which explicitly pre- 
serves the positive semidefinitness of the density matrix 
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It follows from the condition Tr pout — 1 that 

i 

where / denotes the identity operator. Further we can 
expand At in some complete operator basis Aj, 



Ai — 7 , CijAj 



(5) 



If we deal with N level system |i), i = 0, . . . , iV — 1, then 
it is natural to choose the N"^ basis operators as 
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but other constructions are possible. Inserting (0) into 
(0), we find that 



From these iV^ linear constraints one can easily express 
N'^ real parameters in terms of the remaining N^ — N'^ 
ones and thus achieve minimal parameterization. 

The relation between x and Q can be found by com- 
paring Eqs. d) and (g). 
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This formula simplifies considerably if the basis (|6|) is 
chosen Q^^ — XiN+k,jN+i- To provide an explicit exam- 
ple, let us consider a two- level system (single qubit). The 
matrix x can be expressed in terms of ^f- as follows. 
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Gf^. Thus X is 



and the constraints (H) yield G^\ = ^ki 

parameterized by 16 — 4 = 12 real parameters that can 

be collected in a vector 



G = (gO° , ^o^o' , Re goi , Im Cjoi , Re tjo? , Im G^"^ , 

Re gi? , Im G}^ , Re ^0°!^ , Im G^k > Re ^o'l , Im G}k ) (13) 

Note that G\] = (^!i )* since x is hermitian. Additional 
constraints on G\\ follow from the positive semidefinite- 
ness of X- All four main subdeterminants of the matrix 
( |l^ ) should be non-negative. This can be easily checked 
for each G where the likelihood function (0) is evaluated. 
If (O) is not positive semidefinite, then one may simply 



put C[Q] = 0. The maximum of £ can be found for exam- 
ple with the help of downhill-simplex algorithm. In case 
of 2 level system it is sufficient to search for the maximum 
in the finite volume subspace of 12 dimensional space. 

Alternatively, one can find the maximum by setting 
to zero all derivatives of C[G] with respect of Gij ■ It is 
convenient to work with the log-likelihood function. The 
constraints (^ must be incorporated by introducing iV^ 
(complex) Lagrange multipliers Amr; 
arrives at 
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Eqs. (y^ and (^^ represent a system of TV* -t- N"^ non- 
linear equations which must be solved for TV* elements 
Q^^ and N"^ Lagrange multipliers Xmn- On inserting the 
explicit expression for the likelihood function (0) into Eq. 
( |l4| ) one obtains, 






(15) 



where we have introduced 
p^ = TrE A,p|r) Aln('")] = Tr [iv'-^-^ G g'^^) . (16) 

i 

As follows from Eq. (|lj) A is positive definite hermitian 
matrix. 

The extremal equation may be rewritten to the form 
suitable for iterative solution. Multiplying eq. (^ by 
{X~^)inGac ^''^d summing over a,k,l, one gets 
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Convenient form of Lagrange multipliers A,; 
found by inserting eq. (^ into ( pli| ) 

m ^™ a,k,p 



(17) 
may be 

(18) 



The system of nonlinear equations ( |17| ) and (18) for the 
elements of Q can be conveniently solved by repeated 
iterations. 

The theory may be formulated in terms of the opera- 
tors Ai, A\. It is helpful to define a hermitian operator 



A 



^ ATO„|m)(n|. 



(19) 



The maximum of log-likelihood function can be formally 
found as the relation 
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On performing the differentiation with respect to A\ , and 
solving for Ai, we obtain 






(21) 



Next we multiply (|2^) from the left by operator AJ and 
sum over i. Taking into account the constraint (^, we 
find 






(22) 



which is equivalent to (|l^) . Notice that Tr A — ^„^ fm = 
1. 

The procedure of Max-Lik estimation may be inter- 
preted as a generalized measurement. To show this ex- 
plicitly, let us put A: = Hn the relation ( |l5| ) and add all 
the elements over k 
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Since all the traces are equal to 1, this relation reads in 
the operator form 
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This is nothing else as the closure relation for renormal- 
ized positive valued operator measures 
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Moreover, in spite of the fact that the relation used by 
standard reconstructions Pm — fm cannot be fulfilled 
in general, the analogous relation for the renormalized 
POVM is identically true 
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This indicates the privileged role of Max-Lik estimation 
in analogy with the quantum state estimation EG] . Max- 
Lik estimation represents a genuine quantum measure- 
ment. Properties of a quantum black box are determined 
using the closure relation ( p4| ) for a POVM, expectation 
values of which are the registered data (|25|). 

In the rest of the paper we demonstrate the feasibility 
of our approach by means of Monte Carlo simulations for 
two- level system (a single qubit). We shall consider spin 
1/2 system. The detector D shown in Fig. 1 is Stern- 
Gerlach apparatus measuring the spin projections along 
one of three axes x, y, z. We further assume that gin is 
prepared in one of six eigenstates | tj), | \j) of the spin 
projectors (Pauli matrices) CTj, j = a;, j/, z, CTj | \j) = \ tj), 



and 
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Ij). We choose the basis |0) — \ |z) 



and |1) = I 1z). Each of the six input states is prepared 
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FIG. 2. Reconstructed elements of the superoperator Q 
plotted in the form of the vector G. Bars correspond to the 
Majc-Lik estimation (black), linear inversion (grey), and ex- 
act values (hollow). Missing hollow bars indicate the zero true 
values. The superoperator describes the process of damping, 
r,! = 0.5 and V i_ = 0.75, N = 20. 



SA/" times. At the output, one measures M times the spin 
along each of the three axes x, y, z. The corresponding 
six projectors read Uj = \j){j\, j ^U, ix,1y, iy,U, iz ■ 
Let fjk denote the relative frequency of projections to the 
state \k) measured for the input state \j). The likelihood 
function can be expressed as product of 36 terms. 
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where j, k e {U, ix, Ty, iy, h, iz}- 

In our simulations, the black box of the Fig. 1 corre- 
sponds to the damping of /?!„, 
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Here 2r^ > Tji > are transversal and longitudinal 
decay parameters. The elements of reconstructed super- 
operator are depicted in the Fig. 2. The solution was 
obtained by iterations of eqs. ( pTJ ) and ([l^). For the to- 
tal amount of 360 measurements the Max-Lik estimate 
(black) is very close to the exact values Q (hollow). No- 
tice that Max-Lik provides always physically sound re- 
sult, on the contrary to the linear inversion (grey). 

Properties of transforming systems are of interest in 
any physical theory. The developed formalism shows 
how to describe it as a genuine quantum measurement. 
Quantum systems consisting of spins, two entangled or 
three entangled (GHZ) qubits are tractable due to their 
low dimensionality. Proper and full quantum description 
of possible transformations of such systems is, however, 
more advanced, since it is characterized by 12, 240, or 
even 4032 parameters. 
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